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We investigate the finite-size fluctuations due to quenciied disorder of the critical temperature of 
the Sherrington-Kirkpatrick spin glass. In order to accomplish this task, we perform a finite-size 
analysis of the spectrum of the susceptibility matrix obtained via the Plefka expansion. By exploiting 
results from random matrix theory, we obtain that the fluctuations of the critical temperature are 
described by the Tracy- Widom distribution with a non-trivial scaling exponent 2/3. 
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I. INTRODUCTION 

The characterization of phase transitions in terms of 
a non-analytic behavior of thermodynamic functions in 
the infinite-size limit has served as a milestone^"^ in 
the physical understanding of critical phenomena. In 
laboratory and numerical experiments the system size is 
always finite, so that the divergences that would result 
from such a non-analytical behavior are suppressed, 
and are replaced by smooth maxima occurring in the 
observation of physical quantities as a function of the 
temperature. In disordered systems the pseudo-critical 
temperature, defined as the temperature at which this 
maximum occurs, is a fiuctuating quantity depending 
on the realization of the disorder. A question naturally 
arises: can the fluctuations of the pseudo-critical tem- 
perature be understood and determined with tools of 
probability theory? Several efforts have been made to 
study the fluctuations of the pseudo-critical temperature 
for disordered finite-dimensional systems^^^ and their 
physical implications. For instance, recently Sarlat 
et &\}^ showed that the theory of finite-size scaling, 
which is valid for pure systems, fails in a fully-connected 
disordered models because of strong sample-to-sample 
fluctuations of the critical temperature. 

The Extreme Value Statistics of independent random 
variables is a well-established problem with a long 
history dating from the original work of Gumbel^^. while 
less results are known in the case where the random 
variables are correlated. The eigenvalues of a Gaussian 
random matrix are an example of strongly-correlated 
random variables^^. Only recently, Tracy and Widom 
calculated^'^"^^ exactly the probability distribution of 
the typical fluctuations of the largest eigenvalue of 
a Gaussian random matrix around its mean value. 
This distribution, known as Tracy- Widom distribution, 
appears in many different models of statistical physics, 
such as directed polymers^^'^^ or polynuclear growth 
models^^, showing profound links between such different 
systems. Conversely, to our knowledge no evident 
connections between the Tracy- Widom distribution and 
the physics of spin glasses have been found heretofore'^". 



The purpose of this work is to try to fill this gap. We 
consider a mean-field spin glass model, the Sherrington- 
Kirkpatrick (SK) model^^, and propose a definition of 
finite-size critical temperature inspired by a previous 
analysis^. We investigate the finite-size fluctuations of 
this pseudo-critical temperature in the framework of Ex- 
treme Value Statistics and show that the Tracy- Widom 
distribution naturally arises in the description of such 
fluctuations. 



II. THE MODEL 

The SK model^^ is defined by the Hamiltonian 

J N N 

H[{Si},{x,j}] = -—— x,jS,Sj + 'Yh,Si (1) 



where Si = ±1, the couplings {xij}i:^j=i^... = 
{x}, Xji = Xij\/i > i are distributed according to nor- 
mal distribution with zero mean and unit variance 



P{x) 
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and J is a parameter tuning the strength of the interac- 
tion energy between spins. 

The low-temperature features of the SK model have 
been widely investigated in the past and are encoded 
in Parisi's solution^^"^^, showing that the SK has a 
finite-temperature spin glass transition at Tc = J in the 
thermodynamic limit N ^ oo. The critical value Tc can 
be physically thought as the value of the temperature 
where ergodicity breaking occurs and the spin glass 
susceptibility diverges^^"^^. 

While Parisi's solution has been derived within the 
replica method framework, an alternative approach to 
study the SK model had been previously proposed by 
Thouless, Anderson and Palmer (TAP)'^^. Within this 
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approach, the system is described in terms of a free- 
energy at fixed local magnetization, and the physical fea- 
tures derived in terms of the resulting free-energy land- 
scape. Later on, Plefka^^ showed that the TAP free- 
energy can be obtained as the result of a systematic ex- 
pansion in powers of the parameter 



a = 



iVl/2' 



where /3 is the inverse temperature of the model. This a- 
expansion, known as Plefka expansion, has thus served 
as a method for deriving TAP free energy for several 
class of models, and has been extensively used in several 
different contexts in physics, from classical disordered 
systems^""'^^ , to general quantum systems^^"'^^ . It is a 
general fact that, if the model is defined on a complete 
graph, the Plefka expansion truncates to a finite order in 
a, because higher-order terms should vanish in the ther- 
modynamic limit. In particular, for the SK model the 
orders of the expansion larger than three are believed''^ 
to vanish in the limit TV — > oo, in such a way that the 
expansion truncates, and one is left with the first three 
orders of the a-series, which read 
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where = {Si) is the local magnetization, i. e. the 
thermal average () of the spin Si performed with the 
Boltzmann weight given by Eq. (1) at fixed disorder {x}. 

In the thermodynamic limit N ^ oo, for temperatures 
T > Tf. the only minimum of /?/({m}, l3) is the paramag- 
netic one TOj = Vi. Below the critical temperature, the 
TAP free energy has exponentially-many different min- 
ima: the system is in a glassy phase. In this framework, 
the phase transition at Tc can be characterized by the 
inverse susceptibility matrix, which is also the Hessian of 
/ 
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The inverse susceptibility matrix in the paramagnetic 
minimum at leading order in A'^ is: 



(5) 



Random-matrix theory states that the average density of 
eigenvalues of x 
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has a semi-circular shape^^ on a finite support 
[—2\/N,2\/N], where E^. denotes expectation value 
with respect to the random bonds {x}, and Ai({x}) 
is the i-th eigenvalue of x. Eq. (6) is nothing but 
the density of eigenvalues of the Gaussian Orthogonal 
Ensemble (GOE) of Gaussian random matrices^^'-^^. 

Due to self-averaging properties, the minimal 
eigenvalue of fix~^ in the paramagnetic minimum is 
A = (1 - pjf. This shows that, for T > \ is strictly 
positive and vanishes at Tc, implying the divergence^^ 
of the spin glass susceptibility l//3^Tr[x^]. Since A is 
also the minimal eigenvalue of the Hessian matrix of /?/ 
in the paramagnetic minimum, we deduce that this is 
stable for T > Tc and becomes marginally stable at Tc- 

This analysis sheds some light on the nature of the 
spin glass transition of the SK model in terms of the 
minimal eigenvalue A of the inverse susceptibility matrix 
(Hessian matrix) in the thermodynamic limit. In this 
paper we are intended to generalize such analysis to 
finite-sizes, where no diverging susceptibility neither 
uniquely-defined critical temperature exist, and the min- 
imal eigenvalue A acquires fluctuations due to quenched 
disorder. We show that a finite-size pseudo-critical 
temperature can be suitably defined and investigate its 
finite-size fiuctuations with respect to disorder. As a 
result of this work, these fluctuations are found to be 
described by the Tracy- Widom distribution. 

The rest of the paper is structured as follows. In Sec- 
tion HI, we generalize Eq. (5) to finite sizes, in the 
simplifying assumption that the Plefka expansion can be 
truncated up to order , which is known as the TAP ap- 
proach. We then study the finite-size fluctuations of the 
minimal eigenvalue A of the susceptibility matrix, and 
show that they are governed by the TW distribution. 

In Section IV, we extend this simplified approach by 
taking into account the full Plefka expansion, by perform- 
ing an infinite re-summation of the series. 

Hence, in Section V, we give a suitable definition of 
a finite-size pseudo-critical temperature, and show that 
its fluctuations are governed by the TW distribution. In 
Section VL this result is discussed in the perspective of 
generalizing it to more realistic spin glass models. 



III. FINITE-SIZE ANALYSIS OF THE 
SUSCEPTIBILITY IN THE TAP 
APPROXIMATION 

In this Section we study the finite-size fluctuations 
due to disorder of the minimal eigenvalue of the inverse 
susceptibility matrix /3x~^ ^■t the paramagnetic mini- 
mum rrii = OVi, by considering the free energy / in the 
TAP approximation, Eq. (3). We want to stress the fact 
that large deviations of thermodynamics quantities of 
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the SK model have been already studied heretofore. For 
example, Parisi et al. have studied^'' '^""^ the probability 
distribution of large deviations of the free energy within 
the replica approach. The same authors studied the 
probability of positive large deviations of the free energy 
per spin in general mean-field spin-glass models'*^, and 
showed that such fluctuations can be interpreted in 
terms of the fluctuations of the largest eigenvalue of 
Gaussian matrices, in analogy with the lines followed in 
the present work. 

Back to the TAP equations (3), the inverse suscepti- 
bility matrix in the paramagnetic minimum for finite TV 
reads: 



-ax,, + <5„- (1 + + (7) 

V A* 



where 



According to Eq. (8) , is given by the sum of Af — 1 
independent identically-distributed random variables xfj . 
By the central limit theorem, at leading order in N the 
variable is distributed according to a Gaussian distri- 
bution with zero mean and variance 2 
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where Pn{z\ = z) denotes the probability that z\ 
equal to z at finite size A^. 



We set 



S,j (1 + + aM,. 



(10) 



According to Eq. (8), the diagonal elements of My 
are random variables correlated to out-of-diagonal 
elements. The statistical properties of the spectrum 
of a random matrix whose entries are correlated to 
each other has been studied heretofore only in some 
cases. For instance, Staring et al.''^^ studied the mean 
eigenvalue density for matrices with a constraint im- 
plying that the row sum of matrix elements should 
vanish, and other correlated cases have been investigated 
both from a physical"*** and mathematical^^ point of view. 

In recent years, a huge amount of results has been ob- 
tained on the distribution of the minimal eigenvalue of a 
N X N random matrix drawn from Gaussian ensembles, 
such as GOE. In particular, Tracy and Widom*^"*^ de- 
duced that for large A^, small fluctuations of the minimal 



eigenvalue Agoe of a GOE matrix around its leading- 
order value —2\/N are given by 
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GOE 
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Ai/6 



0GOE, 



where 0goe is a random variable distributed according 
to the Tracy- Widom (TW) distribution for the GOE en- 
semble PGOE (</>)• It follows that for j3J = 1 if ^2 '^^^ 
independent on {z}, the matrix Mij would belong to the 
GOE ensemble, and the minimal eigenvalue A of /3x~^ 
would define a variable according to 



A = 



1 



Ar2/3' 



(12) 



and (f) would be distributed according to the TW 
distribution PGOE (</>)• 

As shown in Appendix A, this is indeed the case for 
z\, which can be treated, at leading order in A^, as a 
random variable independent on Xij. The general idea 
is that Z2 is given by the sum of A^ — 1 terms all of the 
same order of magnitude, and only one amongst these 
A^ — 1 terms depends on Xij. It follows that at leading 
order in A^, z\ can be considered as independent on Xij. 
(8) Since in Eq. (7) Z2 is multiplied by a sub-leading factor 



1/v A, in Eq. (7) we can consider Z2 at leading order in 
A^, and treat it as independent on Xij. 

To test this independence property, we set /3 J = 1, gen- 
erate numerically 5^1 samples of the N x N matrix 
/3x^*, and compute the average density of eigenvalues of 
/3x^*, defined as in Eq. (6), together with the distri- 
bution of the minimal eigenvalue A for several sizes A^. 
The eigenvalue distribution pnW as a function of A is 
depicted in Fig. 1, and tends to the Wigner semicircle 
as A^ is increased, showing that the minimal eigenvalue 
A tends to as A^ — )- 00. 

The finite-size fluctuations of A around are then in- 
vestigated in Fig. 2. Defining (p in terms of A by Eq. 
(12), in Fig. 2 we depict the distribution pn{4>) of the 
variable </> for several sizes A^, and show that for increas- 
ing A^, Pn{,4') approaches the TW distribution pgoe (</>)• 
Let us introduce the central moments 

/if = IEjv [</>], 

/if = EJv[(<^-IE^^[(^])*]V^> 1 
of Pat (0), and the central moments 

/ifOE^EcOE [</>], 

/if = Egoe[(</' - IEgoeM)1 V^ > 1 
of the TW distribution, where 



EatH = j d(bpNW; 
1EgOe[-] = / d(t)pGOKi<P) 
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FIG. 1: Density of eigenvalues p]v(A) of the matrix j3x ^ 
for N = 64, 128, 256, 512, 1024 (in red, blue, yellow, green, 
violet respectively), j3J = 1 and S = 16 x 10'^, and Wigner 
semicircular law /9sc(A) — l/(27r)-\/4 — (2 — A)^ (black) as a 
function of A. Pn{^) approaches psc(A) as is increased. 



In the inset of Fig. 2 we depict for several sizes A'' 
and iif^^ as a function of i, showing that /if' converges 
to /if as TV is increased. 

In Figure 3 this convergence is clarified by depicting 
A/if = (/if - /if °^)/Aip°'^ for several values of i > 1 as 
a function of N . Afif is found to converge to for large 
A''. In the inset of Fig. 3 we depict A/if as a function 
of N, showing that the convergence of the first central 
moment with is much slower than that of the other 
central moments. It is interesting to observe that a 
slowly-converging first moment has been recently found 
also in experimental^^ and numerical''^^ data of models 
of growing interfaces where the TW distribution appears. 

The analytical argument proving the independence 
property of has been thus confirmed by this numerical 
calculation. Hence, the main result of this Section is that 
the finite-size fluctuations of the minimal eigenvalue of 
the susceptibility matrix /3x~^ in the TAP approxima- 
tion for pj = 1 are of order N^"^/^ and are distributed 
according to the TW law. These fluctuations have 
already been found to be of order N~'^/^ in a previous 
work^^, and more recently reconsidered''^^ , following an 
independent derivation based on scaling arguments, even 
though the distribution has not been worked out. Our 
approach sheds some light on the nature of the scaling 
A^-2/3, which is non-trivial, since it comes from the 
7V~i/^-scaling of the TW distribution, which is found 



to govern the fluctuations of A. Moreover, the fact that 
we find the same scaling as those found in such previous 
works can be considered as a consistency test of our 
calculation. 

We now recall that both the derivation of this Section 
and the previously-developed analysis of Bray and 
Moore^* rely on the TAP approximation, i. e. neglect 
the terms of the Plefka expansion (13) of order larger 
than 2 in a. As we will show in the following Section, 
these terms give a non-negligible contribution to the 
finite-size corrections of the TAP equations, and so to 
the finite-size fluctuations of the critical temperature, 
and thus must be definitely taken into account in a 
complete treatment. 



IV. FINITE-SIZE ANALYSIS OF THE 
SUSCEPTIBILITY WITHIN THE FULL PLEFKA 
EXPANSION 

In this Section we compute the inverse susceptibility 
matrix /3x~^ by taking into account all the terms of the 
Plefka expansion, in the effort to go beyond the TAP ap- 
proximation of Section III. Notwithstanding its apparent 
difficulty, here we show that this task can be pursued by 
a direct inspection of the terms of the expansion. Indeed, 
let us formally write the free-energy / a a series^^ in a, 

oo 

/({m},/3) = ^a"/„({m},/3). (13) 

For n < 3, the /„s are given by Eq. (3). For n > 3, /„ 
is given by the sum of several different addends"^^, which 
proliferate for increasing n. It is easy to show that at 
leading order in A^, there is just one term contributing to 
/„, and that such term can be written explicitly as 

il>--->i„-l 

x{l~ml) X ••• X (l-m^^ J. 

It follows that by plugging Eq. (14) in Eq. (13) and com- 
puting /3x~^ for rrii = 0, one obtains a simple expression 
for the inverse susceptibility at the paramagnetic solution 



n—3 ii>--->in^ 



-ax^ + (l + /SV^) + 8, 



(15) 
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FIG. 2: Distribution pnW for = 128, 256, 512, 1024, 2048, 4096 (in red, yellow, blue, brown, violet, green respectively) and 
10^ < S < 4 X 10^ samples, and the Tracy- Widom distribution pgoe(</)) for the GOE ensemble (black), as a function of (j)- For 
increasing A'^, PN{<t>) approaches Pgoe('^), confirming the asymptotic independence of the diagonal elements (11) by each of 
the off-diagonal elements Xij for large N. Inset: nf for sizes A'' = 128, 256, 512, 1024, 2048, 4096 (in red, yellow, blue, brown, 
violet, green respectively), 10^ < 5* < 4 x 10^, and nf^^ (black) as a function of i > 1. 



where 



z„ = z — ■-, 



N- 



Xii^Xi^i^ ■ ■ ■ Xi^_^i, Vn > 2. (16) 



il>--->i„-l 



According to Eq. (16), one has that at leading order in 

N 



E,[z;j =0Vn>2, 
E,[(z;)2] = lVn>2, 



(17) 



where in the second hne of Eq. (17) the multiple sum 
defining has been evaluated at leading order in A^. 

We observe that the random variables and Xjk in 
Eq. (15) are not independent, since each depends on 
the bond variables {x}. Following an argument similar 
to that given in Section III for z|, we observe that, by 
Eq. (16) and at leading order in A^, is given by a sum 
of 0{N'"-^^) terms which are all of the same order of 
magnitude. Each term is given by the product of n — 1 
bond variables xu-^Xi-^i^ ■ ■ ■ Xi^_-^i forming a loop passing 
by site i. For any fixed fc and n, only 0(A^"~^) 
terms amongst the 0{N'^~^) terms of z^ are entangled 
with the random bond variable Xjk- It follows that at 
leading order in iV, z^ can be considered as independent 
by Xjk- Since the sum in the second line of Eq. (15) 
has a \/^/N factor multiplying each of the z*jS, we can 
consider the z^ at leading order in TV. Hence, in Eq. 
(15) we can consider each of zj^s as independent on Xj^. 



In Appendix B we show that at leading order in N 
the distribution of z^ is a Gaussian witlr zero mean and 
unit variance for every i and n > 2, while in Appendix C 
we show that at leading order in A^ tire variables {z'^}n,i 
are mutually independent. Both these predictions are 
confirmed by numerical tests, illustrated in Appendix B 
and C respectively. 

Hence, at leading order in A'^ the term in square brack- 
ets in Eq. (15) is nothing but the sum of independent 
Gaussian variables, and is thus equal to a random vari- 
able cr X , wlrere Q is Gaussian with zero mean and unit 
variance, and 



2(/3J)^ 



(n- 1)! 



= 2(/3J)2{2(e('3'^'' - 1) - (/?J)2} 
It follows that Eq. (15) becomes 



axij + 6ij [ 1 + /?- J" + ^=Cj 



N 



where 



a 



(18) 



(19) 



Because of the additional diagonal term in Eq. (19), 
the matrix x'.^j does not belong to the GOE ensem- 
ble. Notwithstanding this fact, it has been shown by 
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FIG. 3: Relative difference A/if between the central mo- 
ments fif of the distribution pjv(<^) for lO'' < 5* < 4 x 10^ 
and the central moments nf^^ of the Tracy- Widom distribu- 
tion as a function of iV = 128, 256, 512, 1024, 2048, 4096, for 
i = 2,3,4,5 (in red, blue, black, orange respectively). For 
increasing A*', fif approaches fif'"'^, confirming the asymp- 
totic independence of Z2 by each of the off-diagonal elements 
Xij for large N. Inset: relative difference of the first central 
moment A/if as a function of A'' (brown). A/if approaches 
very slowly as A'' is increased. 



Soshnikov'^*' that the presence of the diagonal elements 
in Eq. (19) does not alter the universal distribution of 
the maximal eigenvalue of x'^j, which is still distributed 
according to the TW law. Hence, denoting by A the min- 
imal eigenvalue of Px^^, we have 



X = {l-I3jf + 



0GOE, 



(20) 



where (/)goe is a random variable depending on the 
sample Xij , and distributed according to the TW law. 



In this Section we have calculated the inverse sus- 
ceptibility matrix by considering the full Plefka 
expansion. In this framework additional diagonal terms 
are generated that were not present in the TAP approx- 
imation. These additional terms can be handled via a 
resummation to all orders in the Plefka expansion. As 
a result, we obtain that the fluctuations of the minimal 
eigenvalue A of the susceptibility f3x~^ ^re still governed 
by the TW law, as in the TAP case treated in Section III. 



inverse susceptibility matrix, i. e. the Hessian matrix 
of /3/({m},/3) evaluated in the paramagnetic minimum 
TOi = 0, is a function of the temperature and of a 
quantity (/)goe , which depends on the realization of the 
disorder {x}. Since the TW law, i. e. the distribution 
of ^GOE, has support for both positive and negative 
values of (ficOE, the subleading term in Eq. (20) can 
be positive or negative. Accordingly, for samples {x} 
such that 0GOE < 0; there exists a value of /3J sa 1 
such that X{/3J) = 0, in such a way that the spin-glass 
susceptibility in the paramagnetic minimum diverges. 
This fact is physically meaningless, since there cannot 
be divergences in physical quantities at finite size. This 
apparent contradiction can be easily understood by 
observing that if A(/3J) = 0, the true physical suscep- 
tibility is no more the paramagnetic one, but must be 
evaluated in the low-lying non-paramagnetic minima 
of the free-energy, whose appearance is driven by the 
emergent instability of the paramagnetic minimum. 

According to this discussion, in the following we will 
consider only samples {x} such that 0goe > 0. For these 
samples, the spectrum of the Hessian matrix at the para- 
magnetic minimum has positive support for every tem- 
perature: the paramagnetic solution is always stable and 
the paramagnetic susceptibility matrix x is physical and 
finite. We define a pseudo-inverse critical temperature 
/3cJ as the value of (3 J such that A has a minimum at 



dX 



dpj 



-2(1 - /3, J) + ^0GOE (21) 



where in the second line of Eq. (21), Eq. (20) has been 
used. This definition of pseudo-critical temperature has 
a clear physical interpretation: the stability of the para- 
magnetic minimum, which is encoded into the spectrum 
of the Hessian matrix /3x~^, has a minimum at /3 = /3c. 
According to Eq. (21), the finite-size critical temperature 
/3c is given by 



/3c J = 1 - 



1/2 

iV2/3 



(22) 



where 0goe depends on the sample {x}, and is dis- 
tributed according to the TW law. 



V. FINITE SIZE FLUCTUATIONS OF THE 
CRITICAL TEMPERATURE 

We can now define a finite-size critical temperature, 
and investigate its finite-size fluctuations due to disorder. 

In the previous Sections we have shown that for a 
large but finite size the minimal eigenvalue of the 



Eq. (22) shows that the pseudo-critical temperature 
of the SK model is a random variable depending on 
the realization of the quenched disorder: finite-size 
fluctuations of the pseudo-critical temperature are of 
order TV"^/^, and are distributed according to the TW 
law. This has to be considered the main result of this 
paper. 
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VI. DISCUSSION AND CONCLUSIONS 

In this paper, the finite-size fluctuations of the critical 
temperature of the Sherrington-Kirkpatrick spin glass 
model have been investigated. The analysis is carried 
on within the framework of the Plefka expansion for 
the free-energy at fixed local magnetization. A direct 
investigation of the expansion shows that an infinite 
resummation of the series is required to describe the 
finite-size fiuctuations of the critical temperature. By 
observing that the terms in the expansion can be treated 
as independent random variables, one can suitably 
define a finite-size critical temperature. Such a critical 
temperature has a unique value in the infinite-size limit, 
while exhibits fiuctuations due to quenched disorder 
at finite sizes. These fiuctuations with respect to the 
infinite-size value have been analyzed, and have been 
found to be of order iV~^/^, where N is the system size, 
and to be distributed according to the Tracy- Widom 
distribution. 

An analogous role of the TW distribution in the de- 
scription of the critical properties of a physical system 
has also recently been clarified by Forrester et al.^^, 
showing that the TW law describes the finite-size 
corrections of the free-energy of a Yang-Mills theory in 
the neighborhood of its critical point. 

The exponent 2/3 describing the fluctuations of the 
pseudo-critical temperature stems from the fact that the 
finite-size fluctuations of the minimal eigenvalue A of the 
inverse susceptibility matrix are of order N~^/^. Such a 
scaling for A at the critical temperature had already been 
obtained in a previous work''®, where it was derived by a 
completely independent method, by taking into account 
only the first three terms of the Plefka expansion. The 
present work shows that a more careful treatment, 
including an infinite resummation of the expansion, 
is needed to handle finite-size effects. The exponent 
2/3 derived by Bray and Moore'** is here rederived by 
establishing a connection with recently-developed results 
in random matrix theory, showing that the scaling 
N~^^^ comes from the scaling of the Tracy-Widom 
distribution, which was still unknown when the paper 
by Bray and Moore"*® had been written. 

As a possible development of the present work, it 
would be interesting to study the fiuctuations of the 
critical temperature for a SK model where the couplings 
are distributed according to a power-law. Indeed, in a 
recent work^^ the distribution of the largest eigenvalue A 
of a random matrix M whose entries Mij are power-law 
distributed as p{Mij) M~^~'^ has been studied. The 
authors show that if > 4 the fluctuations of A are 
of order N^^^^ and are given by the TW distribution, 
while if /i < 4 the fluctuations are of order 7V~^/^"*/^ 
and are governed by Frechet's statistics. This result 
could be directly applied to a SK model with power-law 
distributed couplings. In particular, it would be inter- 
esting to see if there exists a threshold in the exponent 



fi separating two different regimes of the fluctuations of 
Tc. 

Another interesting perspective would be to generalize 
the present approach to realistic spin glass models 
with finite-range interactions. For instance, a huge 
amount of results has been quite recently obtained for 
the three-dimensional Ising spin glass^^"^°, and for the 
short-range p-spin glass model in three dimensions^*, 
yielding evidence for a finite-temperature phase tran- 
sition. It would be interesting to try to generalize the 
present work to that systems, and compare the resulting 
fluctuations of the critical temperature with sample-to- 
sample fluctuations observed in these numerical works. 
Accordingly, the finite-size fluctuations deriving from 
the generalization of this work to the three-dimensional 
Ising spin glass could be hopefully compared with 
those observed in experimental spin glasses^^, such as 
Feo.sMno.sTiOa. 

Finally, a recent numerical analysis^^ inspired by the 
present work has investigated the sample-to-sample fluc- 
tuations of a given pseudo-critical temperature for the SK 
model, which is different from that defined in this work. 
Even though the relatively small number of samples did 
not allow for a precise determination of the probability 
distribution of that pseudo-critical point, the analysis 
yields a scaling exponent equal to 1/3, which is differ- 
ent from that of the pseudo-critical temperature defined 
here. As a consequence, the general scaling features of 
the pseudo-critical temperature seem to depend on the 
actual definition of the pseudo-critical point itself, even 
though different definitions of the pseudo-critical temper- 
ature must all converge to the infinite-size pseudo-critical 
temperature as the system size tends to infinity. As a fu- 
ture perspective, it would be interesting to investigate 
which amongst the features of the pseudo-critical point 
are definition- independent, if any. 
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Appendix A: Proof of the asymptotic independence 

of Xij and Z2 

Here we show that at leading order in N the variables 
Xij and ^2 are independent, i. e. that at leading order in 

N 

Pn{xij = x,z\ = z) ^ Pwixij =x) X pn{z2 = z)- (Al) 
Let us explicitly write the left-hand size of Eq. (Al) as 
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xi, - 1 



(A2) 




FIG. 4: Probability distribution pjv(z| = z) for S — 10^ and 
different values of TV = 64, 128, 256, 512, 1024 (in red, blue, 
yellow, green, violet respectively) together with a Gaussian 
distribution l/\/27re~^ with zero mean and unit variance 
(black), as a function of z. As A*' is increased pn{zI = z) 

converges to l/\/2ne~^ ^"^ , as predicted by the analytical cal- 
culation, Eq. (B4). Inset: zoom of the above plot explicitly 

showing the convergence of pn{z\ — z) to l/\/27re~^ ^'^ as A'' 
is increased. 



where 



denotes the expectation value with 



respect to the probability distributions of the variables 



5 denotes the Dirac delta function, and 



By plugging Eq. (A4) into Eq. (A2) and using Eq. (A3), 
one has 



1 



X. Zn 



X 



xe ■» 

Pix)^e-- 



= pn{xij = x) X pn[zI = z), (A5) 

where in the first line Eq. (A5) we explicitly wrote the 
expectation value with respect to Xij in terms of the prob- 
ability distribution (2), while in the third line proceeded 
at leading order in N , and used Eq. (9). 



Poo{23 = z) X Poc(zi = z') 
P1024(Z3 = Z,zi = z') 




0.2 
0.15 
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0.05 



12- 

N 



(A3) 



Proceeding systematically at leading order in N, the 
second expectation value in the second line of Eq. 
(A2) is nothing but the probability that the variable 

^^^iw '^k^i k^j ^ki ^ 1) is equal to Z2 ■ We observe that 
according to the central limit theorem, at leading order 
in N this probability is given by 



E 



{Xik}k 



1 (!hT 

-^e ■» 



N 



'^ki 



k^i.k^j 



(A4) 



FIG. 5: ^1024(23 — z, zl = z') for S — 10^ samples (red), and 
the A'' — >■ 00-limit of the right-hand side of Eq. (C3) (black), 
as a function of z, z . For large A'^, pn{z\ — z,zi = z') equals 
Pn{zI = z) X pjv(z4 = z'), as predicted by Eq. (C3). Hence, 
at leading order in N the variables Z3 and zl are independent. 



Appendix B: Computation of the probability 
distribution of zl 



Here we compute the probability distribution of 
at leading order in N . Let us define a super index 
L = {ii, . . . , i„_i}, where L stands for 'loop', since L 
represents a loop passing by the site i. Let us also set 
Xl = Xii^Xi^i^ ■■ -Xi^^^i. By Eq. (16) one has 



D! 



^AL,Vn > 2. 



(Bl) 
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We observe that the probabihty distribution of Xl is 
the same for every L. Hence, according to Eq. (Bl), 
is given by the sum of equally distributed random vari- 
ables. Now pick two of these variables, Xl,Xli. For 
some choices oi L,L', Xl and X^r are not independent, 
since they can depend on the same bond variables Xij . If 
one picks one variable Xl, the number of variables ap- 
pearing in the sum (Bl) which are dependent on Xl are 
those having at least one common edge with the edges of 
X^. The number of these variables, at leading order in 
N, is 0(7V"~^), since they are obtained by fixing one of 
the n — 1 indexes ii, - ■ ■ ,in-i- The latter statement is 
equivalent to saying that if one picks at random two vari- 
ables Xl,Xl', the probability that they are correlated is 



Appendix C: Independence of the z^s at leading 
order in 



Let us consider two distinct variables , , and 
proceed at leading order in N. 

Following the notation of Appendix B, we write Eq. 
(16) as 



N- 



0(iV"-^/iV"-i) = 0{N-'). 



(B2) 



(CI) 
(C2) 



Hence, at leading order in N we can treat the ensemble 
of the variables {Xl}l as independent. According to the 
central limit theorem, at leading order in N the variable 



N— 



(n-l)! L 



Xr 



is distributed according to a Gaussian distribution with 
mean EaifX^l =0 and variance 



1 



(ri-1)! 



(n-l)! 



(B3) 



where in Eq. (B3) Eq. (2) has been used. It follows that 
at leading order in N , z"^ is distributed according to a 
Gaussian distribution with zero mean and unit variance 



where L, L' represent a loop of length n, m passing by the 
site i,j respectively. Some of the variables Xl depend 
on some of the variables X^' , because they can depend 
on the same bond variables Xij. Let us pick at random 
one variable Xl appearing in z^, and count the number 
of variables Xl' in z^ that are dependent on Xl- At 
leading order in iV, these are given by the number of 
Xl' having at least one common bond with Xl, and are 
0(iV™~^). Hence, if one picks at random two variables 
Xl,Xl' in Eqs. (CI), (C2) respectively, the probability 
that XltXl' are dependent is 



It follows that z'j and z^ are independent at leading 
order in N, i. e. for N ^ oo 



Z) ^ 



1 



(B4) 



PJV« = ^>4i = = z)xpn{zI^ = z'), (C3) 



where PN{Zn = z) is defined as the probability that z* 
is equal to z at size N. 



where Pn{zI, 



z, z„ 



z') denotes the joint probability 



that z^ equals z and z^^ equals z', at fixed size N. 



Eq. (B4) has been tested numerically for the first few 
values of n: pN{z'n = z) has been computed by generating 
5^1 samples of {a;}, and so of z^. For n = 3, 4, the re- 
sulting probability distribution p^^z^ = z) converges to 
a Gaussian distribution with zero mean and unit variance 
as N is increased, confirming the result (B4). This con- 
vergence is shown in Fig. 4, where pNizl = z) is depicted 
for different values of together with the right-hand side 
of Eq. (B4), as a function of z. 



Eq. (C3) has been tested numerically for n = 3, m = 4: 
Pn{zI = z,z\ = z') has been computed by generating 
a number S* 3> 1 of samples of {x}, and so of z\,z\. 
As a result, the left-hand side of Eq. (C3) converges 
to the right-hand side as N is increased, confirming the 
predictions of the above analytical argument. This is 
shown in Fig. 5, where Piq2a(z\ = z,z\ = z') is depicted 
together with the N — > oo-limit of the right-hand side of 
Eq. (C3) (see Eq. (B4)), as a function of z,z'. 
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